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Abstract 

The  development  of  high-order  methods  (order  of  aeeuraey  >  2nd  order)  on  unstructured 
grids  is  widely  viewed  as  a  major  paeing  item  in  computational  fluid  dynamics  (CFD). 
Efficient  high-order  methods  capable  of  handling  complex  geometries  are  required  to 
compute  vortex  dominated  flows,  and  to  perform  large  eddy  simulation  and  direct 
numerical  simulation  with  complex  configurations,  and  to  predict  aeroacoustic  noise 
generation  and  propagation.  The  primary  objective  of  the  present  research  is  to  develop 
implicit,  multigrid  and  adaptive  solution  algorithms  for  a  promising  high-order  method, 
the  spectral  difference  (SD)  method. 

Several  major  activities  were  carried  out: 

•  The  development  of  a  hp-adaptation  capability  to  capture  both  discontinuous  and 
smooth  flow  features  with  high  efficiency; 

•  Development  of  a  accuracy-preserving  limiter  for  discontinuity-capturing  for  the 
high-order  SD  Navier-Stokes  solver; 

•  Demonstration  of  the  developed  code  for  real  world  unsteady  vortex-dominated  flow 
problems; 

Significant  Achievements 

•  An  hp-grid  adaptation  capability  has  been  developed  for  the  2D  Euler  solver,  and  also 
demonstrated  for  both  steady  and  unsteady  flow  problems.  Shown  in  Figure  1  are  the 
original  and  adapted  grids  and  the  corresponding  density  contours  for  a  Maeh  3  flow 
inside  a  wind  tunnel  with  a  step.  Grid  adaptation  was  used  to  automatically  refine 
regions  of  discontinuity  and  fine  flow  features  to  achieve  high  efficiency. 

•  The  3D  SD  solver  has  been  successfully  used  to  perform  implicit  large  eddy 
simulation  for  flow  over  a  eorrugate  dragonfly  wing,  which  has  nice  lift  performance 
under  low  Reynolds  number  condition.  Preliminary  simulation  results  indicate  the 
flow  solver  is  capable  of  predieting  the  stall  at  high  angle  of  attack.  Shown  in  Figure 
2  is  the  instantaneous  iso-surfaces  of  the  vorticity  magnitude,  and  the  simulation 
agree  well  qualitatively  in  terms  of  flow  topology.  The  computed  lift  coefficient 
history  is  shown  in  Figure  3,  with  comparison  w’ith  the  averaged  experimental  lift 
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coefficient.  Note  that  there  is  very  good  agreement  between  the  present  simulation 
and  experiment,  demonstrating  the  great  potential  of  the  high-order  solver. 

•  A  hierarchical  moment  limiter  is  being  developed  for  high-ordcr  methods.  In  the 
design  of  this  limiter,  we  adopted  a  “troubled  cell”  approach,  in  whieh  eells  requiring 
limiting  arc  detected  first.  The  design  criterion  for  the  deteetor  is  high  efficiency 
bceause  all  the  eells  arc  examined.  We  also  assume  no  detectors  are  perfeet  as  a  laek 
of  resolution  in  a  smooth  region  can  be  easily  mistaken  to  be  a  discontinuity. 
Therefore  we  demand  that  the  limiter  ean  recover  high-aeeuraey  in  ease  a  smooth 
region  is  deteeted.  An  aeeuraey  study  was  performed  using  the  limiter  and  it  was 
confirmed  that  the  limiter  can  aehieve  the  formal  order  of  aeeuraey  for  smooth 
problems.  The  limiter  was  successfully  employed  to  eompute  the  transonic  flow  over 
the  NACA0012  airfoil  as  shown  in  Figure  4. 
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Figure  1.  Density  contours  and  grid  under  adaptive  h-refinement  at  time=4.0; 
(a)  No  adaptation  (base  grid,  26,238  DOFs);  (b)  3  levels  of  adaptation  (190,017 
DOFs  at  last  time  step);  Grid  is  refined  from  base  grid  every  100  time  steps. 
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Figure  2.  Transient  vorticity  distribution  of  dragonfly  wing  at  Angle  of  Attack  of 
16  deg 

3D  Dragonfly  Wing  at  16  Degree  Angle  of  Attack 


Figure  3.  Computed  lift  history  for  dragonfly  wing  at  angle  of  attack  of  16 
degree  with  comparison  to  the  averaged  experimental  measurement 
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Figure  4.  Density  contours  of  transonic  flow  over  the  NACA0012  airfoil,  at  Mach 
=  0.85,  a  =  1  degree,  computed  with  a  hierarchical  moment  limiter 
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